Z bbbb in the light gluino and hght sbottom scenario 
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The hght gluino (12 ~ 16 GeV) and light sbottom (2 ~ 6 GeV) scenario has been used to explain 
the apparent overproduction of 6-quarks at the Tevatron. This scenario also predicts the decay 
Z bbgg where the gluinos subsequently decay into 6-quarks and sbottoms. We show that this 
can contribute to r4i, — T{Z — > bbbb) since most of the sbottoms and 6-quarks arising from gluino 
decay have a small angular separation. We find that while no excess in r4i, is observable due to 
large uncertainties in experimental measurements, the ratio T{Z — > bbgg) /T{Z — > bbbb) can be large 
due to sensitivity to 6-quark mass, the sbottom mixing angle and the gluino mass. We calculate it 
to be in the range 0.05 — 0.41 inclusive of the entire parameter space. 

PACS numbers: 12.60.Jv, 13.87.Ce, 14.65. Fy, 14.80.Ly 



I. INTRODUCTION 

The Standard Model (SM) has been successful in ex- 
plaining a host of experimental observations on elec- 
troweak and QCD phenomenon at LEP as well as hadron 
colliders. But it is generally believed that the SM is an ef- 
fective theory valid at the electroweak energy scale with 
some new physics lying beyond it. Among the leading 
candidates for a theory beyond the SM is the Minimal 
Supersymmetric Standard Model (MSSM) which has 
been extensively studied in the past few decades. 

The MSSM predicts the existence of SUSY partners 
of quarks, gluons and other known particles of the Stan- 
dard Model. These so-called "sparticles" have not been 
observed which has led to speculation that they might be 
too heavy to have observable production rates at present 
collider energies. However, it has been suggested in Q 
that a light sbottom (6i) with a mass of 0(5 GeV) is not 
ruled out by electroweak precision data if its coupling to 
the Z-boson is tuned to be small in the MSSM. Recently 
Berger et al Q have also proposed a light sbottom and 
light gluino (LSLG) model to explain the long-standing 
puzzle of overproduction of &-quarks at the Tevatron jjj . 
In this model, gluinos of mass 12 — 16 GeV are produced 
in pairs in pp collisions and decay almost immediately 
into a sbottom (2 — 6 GeV) and a 6-quark each. The 
sbottom manages to evade direct detection via i?-parity 
violating decays into soft jets of light quarks around the 
cone of the 6-jet. This mechanism is shown to success- 
fully fit the 6-quark transverse momentum distribution 
at NLO level. Alternatively it has been suggested that 
using the latest 6-quark fragmentation functions reduces 
the discrepancy at the Tevatron In this report how- 
ever we will work in the LSLG scenario. 

Recently, there has been some careful re-examination 
of Z-pole precision data in this scenario. QCD correc- 
tions to the Zhh vertex with sbottom and gluino loops 
have been calculated 0. They contribute negatively to 
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Rii and increase in magnitude with the mass of the other 
eigenstate of the sbottom (62). To maintain consistency 
with data, 62 must be lighter than 125 (195) GeV at 2a 
(3<t) level. An extension of that analysis to the entire 
range of electroweak precision data finds that 62 must be 
lighter than 180 GeV at the 5cr level |7|. A 62 in such a 
mass range would have been produced in association with 
a 61 at LEPII energies (upto 209 GeV) via the couplings 

Z6162 and .Z6261. Since such a 62 has not been observed it 
would seem that LEP data disfavors the LSLG scenario. 
However, these constraints can be relaxed because (i) a 
subsequent study of the decay Z — > bibg+bibg has shown 
that it can contribute positively to Rb and therefore 
should overcome at least some of the negative loop ef- 
fects, (ii) a heavier 62 (^ 200 GeV) might be allowed if 
large CP-violating phases are present in the model @ 
and (iii) experimental searches for SUSY particles are 
heavily model-dependent and, to our knowledge, an ex- 
haustive search of LEPII data for a 62 in this particular 
scenario has not been done. 

In addition to ^-precision data, production of 66-pairs 
at the Z-pole via gluon splittin g h as been re-examined 
recently by Cheung and Keung 10]. They calculate the 
contribution of Z — s- qqgg to the process Z — > qqg* bb 
(q — u,d,c,s,b) in the massless q approximation and 
find that the former is only around 4 — 15% of the 
latter. They do not go further and consider the ratio 
T{Z bbgg)/T{Z bbbb) as they expect it to be very 
similar. Like them we note that the process Z bbgg can 
contribute to T{Z — » 6666) via the final states 66666161, 
66666161 and bbbbbibi which have four 6-quarks and two 
sbottoms arising from the gluino decays g — > 661/661. 
This can happen if the 6-quark and sbottom arising from 
gluino decay prefer a small angular separation, so that 
a typical event looks like Z — > 6666. But Z — > bbgg 
can arise not only from the gluon splitting diagrams in 
Fig. nib) but also from "sbottom splitting" diagrams in 
Fig. n^c). We show that the net SUSY process should 
indeed contribute to Z bbbb and the latter diagrams 
significantly enhance the width for this process. We also 
calculate r{Z — > 6666) to leading order in the SM over a 
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range of 6-quark masses. The final result is a wide theo- 
retical range for the ratio r(Z hhgg)/T{Z bbbb) of 
5% to 41% inclusive of the entire parameter space. 



II. CALCULATIONS 
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FIG. 1: Feynman diagrams contributing to (a) Z bbbb and 
{b),(c) Z —> bbgg. Diagrams with gluon/gluino emission off 
the fe-leg and the crossing of identical particles are not shown. 

Feynman rules for the MSSM given by Rosiek jl^ are 
used to evaluate these diagrams. Their formalism allows 
us to write the lighter sbottom mass eigenstate as a su- 
perposition bl = sin 6*6 6l + cos6'6 6_r of the left and 
right-handed states where 6b is the sbottom mixing an- 
gle. This angle appears in the coupling: 



Zbibi cx {-sin^Ob — -sin^ 



'w 



where 6w is the Weinberg angle. However since elec- 
troweak data excludes the process Z ^ b^bi to a high 
precision, the mixing angle must be fine-tuned to make 
the coupling small i.e., — |sin^0vy, |sf,| ~ 0.38 
where the short-hand notation st = sin 6b is used. We 
vary \sb\ in the narrow range 0.30 — 0.45. 

At constant scale we find that the diagrams in Fig. 
Etc) enhance the width for Z bbgg by 10 - 60%. The 



lower limit is obtained for Sb 



-0.30, rrig = 12 GeV 



and the upper for Sb = -1-0.45, rrig = 16 GeV. Default 
values of nib ~ 4.5 GeV, = 4 GeV are used in this 
analysis and variation of these within mf, = 4 — 5.25 GeV 
and m-j^^ —2 — 6 GeV has little effect. We also verify 
that as mjj^ becomes large, the contribution of Fig. ^c) 



therefore choose the invariant mass of the two gluinos, 
rrigg, as the running scale Q since the diagrams in Fig. 
^b) are still dominant. Using a different invariant mass 
such as Q = TTibg only changes Tbg by 2 — 3%. 

In calculating r4f, the 66-pair produced by gluon split- 
ting cannot be isolated due to interference terms between 
crossed diagrams in Fig. ^a), making the off-shellness 
of the virtual gluon indeterminate. This is in contrast 
to the gluon splitting processes Z qqg* bb and 
Z bbg* —>■ qq, q b, where the secondary production 
of 6-quarks does not interfere with primary production 
at leading order Therefore to calculate T^b we first 
find the ratio r4b/r(Z qqg* bb) at constant Q- 
scale, summing the denominator over q = u,d, s, c, b in 
the massless q approximation. Then T^b is evaluated over 
a running Q-scale as follows: 
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The gluon-splitting process has been studied up to 
next-to-leading logarithm (NLL) level in the past and is 
known to be sensitive to the 6-quark mass. Values rang- 
ing from nib = 4.25 GeV ^3 to variation between the 
pole mass and the i?-meson mass i.e. nib = 4.75 — 5.25 
GeV have been used [ll[l|. DELPHI and OPAL [Ulli 
have measured T^b and they use mb — 4.5 — 5.25 GeV in 
their analysis. Clearly, in the absence of a higher order 
calculation there is no choice but to vary rub over a wide 
range. We vary it from the lower limit of the MS value 
rhbirhb) = 4 GeV to the B-meson mass tub = 5.25 GeV 



The canonical strong coupling value as{Mz) = 0.118 
is used as changes to the running of as have been shown 
to be small 0, 0| even in light of the various new effects 
on the hadronic width of the Z in this scenario 0i 0, S 

iimiii. 



Results are given in terms of the ratios: 



i?4f, 



46 



T(Z — !■ hadrons) 



Rbg = 



SM 



T(Z —> hadrons) 



The total hadronic width of the Z is taken to be r{Z — > 
hadrons) = 1.744 GeV dH]. 

III. RESULTS 

Using eqn. we numerically calculate i?4h = (6.06 — 
3.04) X 10""* for 7716 = 4-5.25 GeV. Measurements of the 
same by DELPHI and OPAL yield (6.0 ±1.9 ±1.4) x 10^" 
and (3.6 ± 1.7 ± 2.7) x 10"'' respectively and the uncor- 
related average is given by the Particle Data Group to 
be = (5.2 ± 1.9) X lO"'' 113. Our entire calcu- 

lated range is within 1.2(T of the experimental average. 
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(a) Rb~g vs gluino mass r vs gluino mass 
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(b) Ri,g vs fc-quark mass 
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FIG. 2: Dependence of Rbg on various parameters. If unspec- 



ified, default values of mj = 12 GeV, nib = 4.5 GeV, m 
GeV and st — +0.38 are used. 



bl 



Therefore, the lack of experimental precision does not 
allow us to fix the 6-quark mass any further. The cen- 
tral value of R^^ is obtained for ~ 4.3 GeV vsrhich 
agrees weW with nii, = 4.25 GeV used in to fit the 
full gluon-splitting process. In a similar fashion, Ri,g is 



FIG. 3: r = Rbg/Rib dependence on gluino mass. Upper 
(lower) bound curves are obtained for mt = 5.25 (4) GeV, 



6 (2) GeV and st = +0.45 (-0.30). 



very sensitive to the gluino mass showing a decrease 
by nearly a factor of 3.5 as rrig varies from 12 to 16 GeV 
[Fig. |2Ia)]. On the other hand, variations in 6-quark 
mass and the sbottom mass (2 — 6 GeV) have very little 
effect (~ 5%) [Fig. EJb)]. The effect of varying Sf, within 
the range |sb| = 0.30 — 0.45 is shown in Fig. ^c) for 
nig = 12 GeV. The variation is ^ 30%, increasing with 
gluino mass. Rbg is lower for negative values of Sb but 
increases with |sh| due to constructive interference with 
the gluon-splitting diagrams. Including all parameters, 
we find Rbg = (0.25 - 1.33) x lO'^. 

The total i?4fc + Rbg equals (3.3 - 7.3) x 10"^ for the 
entire parameter space which is still within 1.2cr of the 
experimental value. However the ratio r — Rbg /Rib can 
be quite large, varying from 5 — 41% [Fig. E]. Thus the 
SUSY process can be a significant fraction (4 — 30%) of 
the total events if it cannot be distinguished from the SM 
bbhh decay. 

In this context we now study the structure of Z — s- hhgg 
events. The cumulative final state qqgg has been stud- 
ied in e+e" annihilation at ^/s = 189 GeV fuj], and 
our results are similar. Fig. 0] shows the opening angle 
{cos 6) between final state fe-quarks, gluinos and gluino 
decay products. We see that the prompt 6-quarks tend 
to be well-separated from each other and also from the 
gluinos. We decay a gluino into a 661/661 pair and find 
that the decay products are rather close to each other, 
with cos9 peaking at around 0.8. The two gluinos on 
the other hand appear to have a more or less uniform 
cosine distribution with a slight preference for smaller 
angles. Therefore most events primarily consist of four 
jets containing 6-quarks, with at least three of them well- 
separated, and two of them having sbottoms in or quite 
close to them. If the sbottom hadronizes completely in 
the detector and deposits it's energy as a jet then it might 
be difficult to separate from the accompanying 6-quark 
without a deliberate search. In that case the events would 
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FIG. 4: Opening angles (cos^) between Z-decay products. 
The prompt 6-quarks are simply denoted by & or & while bdec 
and bdec arise from gluino decay. 



look just like Z — > 5565. However, if the sbottom-jet lies 
near the periphery of the 6-jet then it might widen it 
somewhat. A similar situation could arise if the sbot- 
tom undergoes i?-parity violating decays into soft jets: 
6i u + s;c + d;c + s Some of the 5-jets in the 

Rll^ sample might therefore be unusually wide and/or 
have a higher particle multiplicity. However, this can also 
happen due to radiation of hard gluons from 5-quarks in 
Z ^ 6655, and might be difficult to distinguish. The 
remaining small fraction of events where the sbottom(s) 
are well-separated from the 6-jet would compete with the 
background from Z —^ bbbb + j and Z —^ bbbb + 2j events 
and might require careful separation. Finally, if the sbot- 
toms deposit only a fraction of their energy in the detec- 
tor they might still be observed as missing energy. This 
energy is likely to be small though since sbottoms arise 



from gluino decay and tend to be rather soft in compari- 
son to prompt 5-quarks. Events with missing energy can 
also arise from Z bbbb because of subsequent semilep- 
tonic decays B Xlv. While a full detector simulation 
is beyond the scope of this report; given the experimental 
uncertainties and the absence of a deliberate search it is 
possible that sbottoms produced in this scenario escaped 
detection. The SUSY events might therefore lie hidden 
in data for Z bbbb and could be sufficient in number 
that a deliberate search might uncover them. 



IV. CONCLUSIONS 

We show that the light sbottom and light gluino sce- 
nario predicts the SUSY decay Z bbgg. We also show 
that this can contribute to the SM process Z bbbb be- 
cause of gluino decay into pairs of proximate sbottoms 
and 5-quarks. This scenario cannot be constrained by 
an excess in the measured rate due to large experi- 
mental and theoretical uncertainties. However, the ratio 
between the SUSY and SM decays is shown to be signif- 
icant, from 5% to 41% inclusive of the entire parameter 
space. A good fraction of events in the sample might 
therefore arise from the SUSY process and contain sbot- 
tom signatures hidden in and around 5-jets. We suggest 
a model-dependent search for sbottoms in existing LEP 
data in order to constrain this scenario. 
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